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PREFACE 


This  report  is  part  of  RAND's  continuing  interest  in  the  geo* 
physical  effects  produced  by  nuclear  explosions.  Previous  work 
(RM-4225,  RM-4388,  RM-4494,  RM-4858,  RM-4946,  RM-5616)*  has  been  con¬ 
cerned  with  the  geomagnetic  effects  produced  by  such  explosions.  The 
present  study  is  the  start  of  an  investigation  of  some  aspects  of  the 
acoustic  effects.  This  work  was  sponsored  by  the  United  States  Air 
Force  and  the  Advanced  Research  Projects  Agency. 
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SlftMARY 

A  complete  asymptotic  analysis  Is  carried  out  for  the  flow  field 
produced  by  the  Instantaneous  release  of  energy,  at  a  point  on  the 
ground,  In  an  Isothermal  atmosphere.  A  double  Integral  representa¬ 
tion  of  the  flow  Is  constructed  from  Laplace -Hankel  transforms  of 
the  linearized  equations.  An  asymptotic  approximation  to  the 
Integral  Is  obtained  by  two  successive  applications  of  the  method 
of  stationary  phase.  It  Is  found  that  there  are  three  principal 
groups  of  dispersive  waves  behind  the  spherical  acoustic  front.  One 
of  these  groups  Is  contained  In  a  high  frequency  band,  and  the  other 
two  In  a  low  frequency  band.  The  spatial  domain  of  the  low  frequency 
waves  Is  cut  off  by  a  front  (caustic),  the  location  of  which  Is 
calculated.  Some  characteristics  of  the  flow  are  discussed. 
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I.  INTRODUCTION 

The  problem  of  waves  in  a  stratified  atmosphere  has  received 
much  attention  and  many  different  cases  have  been  considered  (cf 
Bibliography)  .  The  case  which  is  analyzed  here  is  one  for  which  a 
complete  asymptotic  analysis  can  be  carried  out  and  an  overall  view 
of  the  resulting  wave  motion  obtained.  The  principal  simplification 
is  that  the  waves  are  small  disturbances  to  an  isothermal  atmosphere. 
The  problem  worked  out  in  detail  is  the  flow  field  produced  by  the 
instantaneous  release  of  energy  at  a  point  in  the  ground  plane  (the 
analysis  is  easily  extended  to  cases  of  energy  release  at  a  point 
above  the  ground)  .  This  case  does  not  seem  to  have  been  discussed 
previously,  although  similar  problems  for  monochromatic  waves  have 
been  discussed. ^ 

The  basic  method  used  for  obtaining  the  asymptotic  approximation 
is  the  repeated  application  of  the  method  of  stationary  phase.  The 
kinematic  aspects  of  the  groups  of  waves  produced  could  have  been 

(2) 

based  on  the  dispersion  relation  and  the  considerations  of  Whitham'  ' 
(3) 

and  Lighthill.  But,  it  turned  out  to  be  more  convenient  for  this 
paper  to  base  all  the  calculations  on  the  Fourier  integrals  derived 
by  Lap  lace  -Hankel  transforms  of  the  linearized  equations.  In  any 
case  there  are  three  principal  groups  of  waves,  one  in  a  high  fre¬ 
quency  band  and  two  in  a  low  frequency  band.  The  spatial  domain  of 

the  low  frequency  waves  is  cut  off  by  a  front  (caustic) .  This 

(4) 

phenomena  was  first  noted  by  Mowbray  and  Rarity'  in  their  study  of 
dispersive  waves  in  an  incompressible  stratified  fluid.  A  similar 
phenomenon  appears  in  the  pattern  of  waves  around  a  moving  ship. 

In  Section  II  the  basic  equations  are  given  and  the  lineariza¬ 
tion  is  carried  out  in  suitable  dimensionless  units.  The  basic 
small  parameter  is  shown  to  be  the  ratio  of  the  energy  released  to 
the  internal  energy  in  a  volume  of  scale  height  dimensions.  In 
addition,  due  to  the  exponential  decay  of  density,  there  is  an  ex¬ 
ponential  growth  of  disturbances  unfavorable  for  the  validity  of  the 
theory.  However,  there  should  still  be  a  considerable  range  in  which 
the  approximation  is  quantitatively  correct  as  well  as  qualitatively 
instructive. 
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In  Section  III  the  application  of  Laplace  transform  in  time  and 
Hankel  transform  with  respect  to  cylindrical  radius  enables  an 
integral  representation  of  the  exact  solution  to  be  obtained.  In 
Section  IV  it  is  shown  how  the  acoustic  spherical  wave  front  is  con¬ 
tained  in  the  representation.  Non-linear  effects  are  probably  most 
important  in  describing  the  amplification  of  the  front  into  a  shock 
wave  and  its  subsequent  propagation.  However,  the  main  residual 
disturbance  is  a  system  of  dispersive  waves  behind  the  front.  By 
studying  the  double  integral  representation,  it  is  seen  that  a 
cylindrical  ground  wave  exists  in  addition  to  propagating  waves  in 
low  and  high  frequency  bands.  Section  V  studies  these  propagating 
dispersive  groups  and  presents  the  final  description  of  the  flow  as 
well  as  some  numerical  results. 
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II.  BASIC  EQUATIONS 

We  consider  motion  in  an  isothermal  atmosphere  above  a  ground 
plane  produced  by  instantaneous  energy  release  at  a  point  on  the 
ground,  the  origin  in  Fig.  2.1.  The  equilibrium  atmosphere  is  thus 
characterized  by  the  usual  exponential  distributions  with  scale 
height  h. 


P0(z) 

* 

0 


(2.1) 


scale  height 


*  *  * 

P  ,  p  ,  T  ■  sea  level  (z 

* 

RT 

r~S  ®  ' 

sound  speed  -  Jy RT  ,  cg  ■ 

less  equations  of  motion, 
gas  are 


-  0)  pressure,  density,  temperature, 

Note  that  c  ■  Jy  c^  where  c  ■  isentropic 
gravity  wave  speed  -  ,/gh.  The  dimension- 
continuity,  momentum,  entropy  for  a  perfect 


+  div  pq  ■  0 


(2.2) 


ff+J-'S-  -Ifc  (2.3, 

<!t  +  h  •  h  «  .  (2.4) 

PT  o’ 

Here  lengths  have  been  referred  to  the  scale  height  h,  velocities  to 

•jg  jdf  fg  <fg 

c  ,  times  to  h/c  ,  pressure  to  P  ,  and  density  to  p  .  The  basic 
parameter  of  the  problem  is 

(v-DQ0 


C 


(2.5) 
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Fig.  2.1  -  Co -ordinate  system,  with  origin  at  tha  point 
of  anargy  ralssse. 
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where  Q  -  energy  released  at  t  ■  0;  e  roughly  measures  the  energy 
o 

release  compared  to  the  internal  energy  stored  in  a  scale  height 
volume.  The  initial  conditions  at  t  ■  0  are 

P  *  p  -  e"z,  q  -  0,  (2.6a) 

4 

the  boundary  conditions  at  z  ■  0  is 

q  -  0.  (2.6b) 

z 

ACOUSTIC  EXPANSION 

The  acoustic  expansion  is  based  on  e  «  1  and  represents  the 
flow  as  small  changes  superimposed  on  the  ambient  state: 

q  -  e  u  +  . ..,  u  -  (u,w)  (2.7) 

where  u  -  radial  component 

w  ■  vertical  component 


P  -  e"z  {l  +  ep  +  ...}  (2.8) 

p  -  e  z  (l  +  9.o  +  ...}  .  (2.9) 


The  equations  of  order  e  that  result  from  (2.2,  2.3,  2.4)  are  thus, 
for  the  case  of  axial  symmetry, 


+  u  + 
t  r 

u  ,  n  /  2  2 

-  +  w  -  w  -  0  r  -  V*  +y 

1  2 

(2.10) 

ut 

-  -Ip 

y  *r 

(2.11) 

wt 

-  -ip 

Y  z  Y 

(2.12) 

Pt  -  Y  at 

+  (y-l)w  -  6(z)6(t)  . 

(2.13) 
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The  initial  and  boundary  conditions  for  the  perturbations  are 

p  -  ct  ■  u  ■  0  at  t  ■  0-  (2.14) 

v  -  0  at  z  -  0,  t>0  (2.15) 

and  the  main  wave  motion  is  assumed  outgoing.  The  characteristic 
surfaces  of  the  system  (2.10,-13)  are  given  by 

*\  -  0,  0*  -  -  0*  -  o  (2.16) 

where  0(r,z,t)  -  const,  on  a  characteristic  surface  in  (r,z,t)  space. 
The  first  factor  in  (2-16)  represents  possible  discontinuity  surfaces 
along  fixed  cylinders  r(z)  in  space  while  the  second  shows  the 
possibility  of  a  spherical  acoustic  front  (the  sound  speed  is  con¬ 
stant  everywhare  in  this  model) . 
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III.  FORMAL  SOLUTIONS 

A  formal  Integral  representation  of  the  solution  is  constructed 
by  a  combination  of  Laplace  and  Hankel  transformation.  Let 


®  _  t  * 

g(s)  -  ]  e~s  g(t)  dt  (Laplace  transform) 

(3.1) 

o 

so  that  the  system  (2.10,13)  has  ^  -•  s,  with  the  zero 

initial 

conditions . 

••  ■  u  —  ■ 

su  +  u  +  —  +  w  -  w  ■  0 

r  r  z 

(3.2) 

■  1  OO 

su  +  —  p  “0 

Y  r 

(3.3) 

^  a  +  sw  +  i  (Pz  -  p)  -  0 

(3.4) 

-  ysa  +  (y-1)w  +  sp  -  ~  6(z). 

(3.5) 

Using  the  Hankel  transforms 

f  (r)  -  J  F(U)>  Jv(uur)u)du),  F(u>)  -  J  f(r)  Jv(Wr)rdrJ  v  >  (-  $ 
o  ° 

(3.6) 


note  that 


-  J  6(w)  6(»)  -  J  J0(w)rdr  -  ^  . 

o  ° 

(3.7) 

intend  to  use  the  usual  inversion  formula 
i  00 

g(t)  -  l/2TTi  J  est  i(s)ds  . 

-i  00 
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Thus ,  let 

00 

(a,w,p)  -  J  (R,W,P)  JQ(kr)  kdk  (3.8) 

o 

oo 

u  =  J  U(k)  (kr)  kdk  .  (3.9) 

o 

The  choices  In  (3,8,9)  are  made  consistent  with  the  overall  symmetry 
of  the  equations  and  insure  that  £i^  ruf  -  0,  so  that  the  behavior 
at  the  axis  of  symmetry  is  regular.  Using  the  properties  of  Jq, 

J^,  etc.  the  basic  system,  (3.2,5)  becomes 


dU 

sR  +  kU  +  ~  -  W  *0 

(3.10) 

■»  -  Tt  '  '  0 

(3.11) 

-  R  +  sW  +  -  (j|£  -  P)  -  0 
y  y  'dz 

(3.12) 

-  y  sR  +  (y-l)W  +  sP  -  ^  6(Z)  . 

(3.13) 

This  system  is  linear  and  second  order  with  respect  to  z,  so  that 
there  are  two  basic  roots  V. ,  \0,  for  solutions  of  the  form 
(RtU,WtP)  ~  e  z.  The  characteristic  determinant  from  (3.10,-13) 
is 

s  k 

0  s 


-ys  0 


X-l 

0 


y-1 


0 

-k/y 

-(\-l) 

Y 


(3.14) 


and  the  characteristic  equation  is 
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(X-l)2  +  (X-l)  -  {s2  +  k2  (1  +  -^J)}  -  0  . 

V  +8 


(3.15) 


The  basic  roots  are 


+  k2(l  + 


v-l 
2  2 
V  s 


(3.16) 


For  definiteness  the  contours  in  the  ®  plane  is  at  first  chosen 
so  that  Re  >0  in  (3.16). 

In  order  to  formulate  "jump"  conditions  it  is  most  convenient 
to  eliminate  (R,U)  from  the  system  (3.10-13)  and  obtain 


-  W  + 


s)P 


(3.17) 


{v2s2  +  (Y-l)}  W  +  ys  (^  -  P)  +  sP  -  ~  6(z)  .  (3.18) 


Integration  of  (3.17)  from  (0-  to  0+)  yields  the  jump  conditions 


[w] 


zK) 


_L.  rp]  -  JL I 

2tty  *  L  J2^)  2tty  S  * 


(3.19) 


Since  the  solution  should  vanish  ahead  of  the  acoustic  wave  front 
dtt  >  t)  only  the  root  (3.16)  need  be  considered  and  we 
write 


X  *  h  -  *i(s,k;Y)  where  i*  +  s2  +  k2(l  +  (3. 


20) 


,2  .  Ill 


2  * 


and 


-10- 


2  2 

Note  that  0  <  P  <  fc;  (y  ■  1.4,  P  ■  .22)  for  realistic  values  of 
y.  Since  W  *  0,  for  z  ■  0-,  the  Jump  condition  (3.19)  provides  the 
boundary  conditions  for  W  as  z  -*  0+  and  the  solution  for  W  is 

W(z;s,k)  ■  j“e^  e~^z  .  (3.21) 

It  follows  from  (3.17),  for  z  >  0,  that 


P(z ;s ;k) 


1_ 

2tt 


2-y 

2v 


+  » 


s2  + 


s  . 


(3.22) 


Applying  the  double  inversion  of  Laplace  and  Hankel  transforms,  we 
obtain  a  formal  representation  of  the  perturbation  pressure  field 


where 


e"^Zp(r,z,t)  “  2n  J  Jo*kr*  PC11*2*1)  kdk 

o 


F(k;z,t) 


JL  P  -8t~**(8»k>z 

2rri  J 

-i® 


^1  + 
2y 

s2  +  k 


s  ds 


(3.23) 


(3.24) 


The  integration  of  (3.24)  is  regarded  as  carried  cut  first,  for  a 
fixed  real  k.  This  integration  over  s  is  carzfed  out,  at  first, 
along  a  path  parallel  to  the  imaginary  axis  to  the  right  of  all 
singularities  in  the  ®  plane,  so  that  the  initial  conditions  are 
satisfied.  The  ©  plane  is  cut  suitably  so  that  Re(^)  >  0  on  the 
contour . 

The  pressure  field  is  the  main  quantity  of  interest  but  all  the 
other  quantities  easily  follow  from  the  system  (3.10-13). 
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IV.  APPROXIMATE  AND  ASYMPTOTIC  EVALUATIONS 

In  this  note  the  main  Interest  Is  centered  on  the  dispersive 
waves  behind  the  acoustic  front.  This  is  the  main  disturbance  which 
contains  eventually  most  of  the  energy  put  in  at  t  *  0. 

However,  first,  a  crude  approximation  valid  near  the  acoustic 
wave  front  comes  from  considering  large  (k,s)  in  (3.23),  to  obtain 


p(r,z,t) 


•  &  1  f  st . 

*  e  — s  e  ds 

4tt  i  -i® 


/  2.2 

®>  -z  Js  +k 

e~JZF~  J°(kr)kdk  • 


(4.1) 


The  second  integral  is  a  special  case  of  an  integral  evaluated  in 
Watson's  Bessel  functions. ^ 


o 


e-z/s2+k2 

J5? 


JQ(kr)kdk 


ZTT 

-s./r  +z 


Jt\zZ 


(4.2) 


so  that  (4.1)  Is 


,  .  e^z  1  f1"  s(t-R) . 

p(r’,'t)  'mm]  **  d“ 


-i® 


(4.3) 


n  t 

where  R  -  *fv  +z  -  spherical  radius.  According  to  the  usual  inter¬ 
pretation  the  Integral  in  (4.3)  is  the  derivative  of  a  delta  function, 
or 


p(r,z,t)  *  2^6,(t“R) 


(4.4) 


6'(t-R)  is  the  acoustic  version  of  the  compression  followed  by  an 
expansion  that  is  felt  near  the  wave  front.  The  intensity  is 
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doubled,  due  to  the  ground  plane,  over  that  of  a  point  source  in  free 

kz 

space.  Further,  the  factor  e  represents  the  amplification  of  the 
waves  proportional  to  pQ(z)  »  according  to  the  usual  local  energy 
cons iderations . 

Asymptotic  integration  of  (3.23,24)  which  is  valid  in  the  main 
wave  zone  for  large  (r,t)  can  be  carried  out  by  repeated  applica¬ 
tions  of  the  method  of  stationary  phase  coupled  with  the  use  of  an 
asymptotic  representation  of  JQ. 

In  order  to  give  definite  meaning  to  the  integral  the  follow¬ 
ing  plan  is  adopted.  The  ©  plane  is  suitably  cut  and  the  integral 
of  (3.24)  for  F(k;z,t)  is  expressed  as  an  integral  along  segments 
of  the  imaginary  axis  (s-iu>) .  That  is,  it  is  expressed  as  a  Fourier 
integral  over  the  frequencies  u>  which,  for  a  given  k,  can  propagate. 
In  doing  this,  the  residues  at  certain  poles  (due  to  presence  of  the 
ground),  have  to  be  evaluated.  Next,  in  the  real  double  Integrals, 
(u>,k) ,  the  order  of  integration  is  changed,  since  the  stationary 
phase  approximation  can  more  easily  be  carried  out  in  the  k-integral. 
However,  this  can  only  be  done  if  the  asymptotic  approximation  for 
JQ  is  used.  After  this  is  completed,  the  answer  is  expressed  as  a 
Fourier  integral.  This  integral  is  again  of  such  a  form  that  the 
method  of  stationary  phase  can  be  used  and  its  qualitative  features 
are  easily  discussed.  (See  5.)  In  addition  some  quantitative  re¬ 
sults  are  also  presented.  Some  details  of  the  procedure  just  out¬ 
lined  are  now  given. 

First  note  that  the  exponent  in  (3.24) 

H(s,k)  -  {s4  +  (k  +  k2)s2  +  02k2  }  (4.5) 

can  be  written  as 

»  -  “  “>*00)  (s 2  +  <4*6> 


so  that  the  integrand  has  branch  points  at  s  -  +  iu>^,  +  iiu^  where 
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*2(k)  “  h  [tt  +  k2  +  20k)  ^  +  (k  +  k2  -  20k)  ^  ]  (4.7) 

^(k)  -  h  [ft  +  k2  +  20k) ^  -  (*  +  k2  -  20k) ^  ]  (4.8) 


u>2  is  a  lower  cut-off  frequency,  u)^  an  upper  cut-off  frequency. 

(See  Fig.  4.1)  In  addition,  the  integrand  of  (3.24)  has  poles  at 
s  =  +  ik,  and  an  essential  singularity  (due  to  p.  in  the  exponent) 
at  the  origin.  It  can  be  verified  that  the  integrand  is  one-valued 
in  the  entire  plane  with  barriers  as  shown,  if  the  principal  branches 
of  the  angles  associated  with  each  branch  point  are  chosen  as  shown 
in  Fig.  4.2.  The  picture  is  drawn  for  a  given  k  and  is  qualitatively 
the  same  for  all  k.  Note  the  arguments  of  p,  at  the  various  points 
on  the  imaginary  axis . 

arg  p,A  -  arg  ^  -  arg  Hg  -  arg  ^  -  tr/2  (4.9) 

arg  -  arg  p>c  -  arg  p^  -  arg  J*G  -  -  tt/2  .  (4.10) 


The  application  of  Cauchy's  theorem  enables  the  integral  of  F  on  the 
original  path  to  be  expressed  as  integrals  on  segments  of  the 
imaginary  axis  as  shown,  plus  (2rri)  x  (sum  of  the  residues)  at  the 
poles  (s  ■  +  ik) .  This  is  true  since  the  integral  on  a  large  circle 
in  the  left  half  plane  vanishes  (Re(s)  <  0).  Also,  the  essential 
singularity  at  the  origin  is  outside  the  contour  used  for  Cauchy's 
theorem,  and  thus  does  not  contribute  a  residue.  Thus 

F(kjz,t)  ■  (  /L  Res)  +  F^  +  F2  +  F-  +  F-  +  F-  +  Fg 

s-  +ik  *  4  J 


(4.11) 
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Fig.  4.2  -  Contour  £or  carrying  out  the  integration  in  the 
8-plane,  with  the  location  of  pole*),  branch 
points,  and  essential  singularity  as  indicated. 
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where 


F^k) 
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.“2  i(0)t+u  z)  -  iU,) 
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and  p  =  ip,^  on  1  »  etc. 


^  =  i  v^C^2"^)  C^2“^2>  =  £/w4-C*+ k2)u)2+p2k2  (4 
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Jar! 


-U) 
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(4 


Noting  that  p,^  =  -  p,^  these  integrals  can  be  combined  to 


lp“  uj  du)  e  2-V  i(u>t-p^z)  2-Y 

f1+F2+f5+F6  *  2^T  f  2  2  t("2y  +  iul)  e 

-oo  (D  -k 


"  (  2y  "  ltii 


"I  >  “>2(k)  . 


(4 
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.14) 


.15) 
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i(u*+U,zs 

t 


.17) 


The  integration  is  carried  out  in  (4.17)  only  over  those  frequencies 
M  -■  (^(k)  which  can  propagate  as  waves.  In  a  similar  way, 
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l(u)t+^z) 


(4.18) 


Thus,  all  these  integrals  can  be  combined  to  a  Fourier  integral  on 
only  those  frequencies  which  can  propagate. 


Fl+F2+F3+F4+F5+F6 


1  r"  duu  f . 
“2niJ  lo»2  — ic.2 1  ** 


2-Y 

2Y 


i(u>t 


•HjZ) 


,2-Y  .  v 


uu|  <  a)L,  |u)|  >  u>2 


(4.19) 


Further,  the  residue  at  s  ■  +  ik  is  easily  calculated: 


(Vi2(ik)  -  *  -  P2  -  (~)2) 

2-Y 

_  .  ... .  2-Y  "  2Y 

Res  (s  -  +  ik)  =■  e 


z  +  ikt 
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(4.20) 


Thus 
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Res 


s"  +ik 


2-Y 

2-Y  "  2Y  Z 

- e 

Y 


cos  kt  . 


(4.21) 


The  pressure  field  is  calculated  by  applying  the  Hankel  transform  to 

(4.19)  and  (4.21)  as  in  (3.23).  The  part  produced  by  the  residue 

at  the  poles  (4,21)  is  called  the  ground  wave  P  and  thus  has  the 

G 

representation 
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iODt+^z) 
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cos  kt  k  dk  . 


(4.22) 


The  Integral  in  (4.22)  is  divergent,  but  can  be  interpreted  as 


sin  kt  JQ(kr)  dk  = 


o 


(t2-r2)3/2 


(4.23) 


Then, 


2-V 

2Y 
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-t 


\(t2-r2)3/2, 


(4.24) 


(4.24)  represents  a  cylindrical  wave,  which  exists  only  behind  the 
spherical  front.  It  is  an  exact  solution  to  the  original  equations 
as  given  in  Lamb^(p.  548).  The  wave  is  excited  from  the  axis 
(r-0)  by  the  vertical  passage  of  the  spherical  acoustic  wave.  The 
remainder  of  the  pressure  field  comes  from  the  Hankel  transform  of 
(4.19).  The  order  of  integration  is  changed  and  the  integral  is  broken 
up  into  low  and  high  frequency  bands  (cf.  Fig.  3.1) 

p  ■  pG  +  Pi  +  PlI  (4.25) 
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Here 


k^o))  -  h  j -j* 


from  (4.14).  These  integrals  can  be  expressed  as  integrals  over 
positive  frequencies  in  low  and  high  frequency  bands. 
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The  asymptotic  formula 
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(4.32) 


is  now  used  in  (4.30,  4.31)  .  This  approximation  is  very  good  for  r 
large  except  for  those  k  near  zero  (in  G^) .  But  for  k  close  to 
zero  the  integrand  is  small  anyway  so  that  the  approximation  of  (4.32) 
is  probably  excellent.  The  apparent  singularity  at  (r-0)  in  (4.32) 
in  fact  disappears  from  the  answer  later.  Thus, 
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Here 

(k)  =  kr  -  Uj(k)z  -  jr  sign  k  0^(-k)  =  -  kr  -  p,^z  +  ^ 

02(-k)  “  -  kr  +  Hj*  +  J 


02(k)  =  kr  +  a^z  -  ^  sign  k 


(4.35) 
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(4.38) 


Now,  the  method  of  stationary  phase  (Ref.  6,  p.  395)  is  applied  to 
the  integrals  in  (4.33)  and  (4.34)  with  the  rapidly  oscillating  ex¬ 
ponents  0.  „(+  k)  .  The  main  contribution  comes  from  wave  numbers 
1  **• 

where  0'(k  )  =0  (main  group).  Details  of  the  calculation  are 
s 

shown  only  for  G^.  The  point  of  stationary  phase  in  0^(k)  is  given 
by 


(4.39) 


dp 

Since  >  0  in  ,  only  0^(k),  (^("k)  =*  -  0^(k)  have  stationary 
points  and  contribute  the  dominant  terms  to  the  answer.  Thus 
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Here  v  is  the  pole  angle  in  physical  space,  r  =  Rsin  cp,  z  ■  Rcos  cp. 
According  to  (4.40),  a  stationary  point  occurs  in  only  for  those 
frequencies  tu  >  8  cos  cp.  According  to  the  method,  the  integrand  is 
evaluated  at  the  stationary  point  except  for  the  exponent  which  is 
approximated  as 

(R-k  )2 

ax(k)  -  —  aj(ks)  (4.4i> 


(subscript  s  denotes  quantities  evaluated  at  the  stationary  point) . 
Due  to  the  rapid  oscillation  away  from  k-kg,  the  limits  of  integra¬ 
tion  can  be  extended  to  »  without  introducing  a  significant  error. 
Note  the  formulas 
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where 
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For  Gjj,  the  only  stationary  points  occur  in  (^(k)  and  O^(-k).  Thus, 
the  stationary  phase  approximation  yields, 
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The  integrals  in  (4.46),  (4.47)  are 
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Using  the  expressions  at  the  stationary  point  tabulated  above,  we 
obtain  for  (4.46,  4.47) 
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(4.49) 

Thus,  the  pressure  in  the  far  field  can  be  written  as  integrals  over 
the  important  frequency  bands  (cf.  4.28,  4.29): 
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(4.50) 
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(4.52) 


The  representation  (4.50,  4.51)  shows  spherical  waves  with  a  spherical 
phase  and  group  velocity  depending  on  angle  cp.  These  integrals  can 
again  be  approximated  by  the  method  of  stationary  phase.  Details 
are  given  in  the  next  section. 


V.  ASYMPTOTIC  PRESSURE  FIELD.  NUMERICAL  RESULTS 


In  this  section,  we  shall  examine  some  of  the  properties  of  the 
pressure  field  given  by  (4.50)  and  (4.51).  The  field  has  been  ex¬ 
pressed  in  terms  of  integrals  of  the  form 

J  f(o),  cos  <p)  du)  (5.1) 

over  two  bands  of  propagating  frequencies,  with  Q(u),cp)  given  by 
(4.52).  For  large  R  and  t,  the  exponentials  are  rapidly  oscillating 
and  the  method  of  stationary  phase  can  be  applied  to  these  integrals. 
The  points  of  stationary  phase,  if  any,  are  given  by 


The  existence  and  number  of  stationary  points  in  each  frequency  band 
is  best  illustrated  graphically. 

We  first  consider  the  high  frequency  band,  Y  <  u>  <  ®.  In 
Fig.  5.1,  we  have  plotted  fi(u>,q>)  vs  u>  for  a  typical  value  of  cp. 

The  function  fi  is  a  monotonically  increasing  function  of  u>,  start¬ 
ing  out  as  ^<d-Y  for  u)  near  Y  and  asymptotically  approaching  u) 
as  -*  «o.  xhe  derivative  ^  decreases  continuously  from  »  to  1 
as  u)  increases  from  Y  to  «».  Since  ^  >  0,  it  is  clear  that  only 
the  exponential  with  the  minus  sign  can  have  a  stationary  point. 

It  is  also  clear  from  the  above  discussion  that  there  will  be  a 
single  stationary  point  in  this  frequency  range  for  all  t/R  such 
that  1  <  t/R  <  «»,  i.e.,  for  all  points  behind  the  spherical  front. 

For  the  low  frequency  band,  8  cos  <  u>  <  8,  the  situation  is 
somewhat  different,  as  can  be  seen  from  Fig.  5.2,  The  function  n 
increases  monotonically  with  uj  in  this  band,  starting  out  as 
«/w-@  cos  for  u)  near  8  cos  <p  and  approaching  *  as 


when 
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3 


<v 


Fig.  5.1  -  High  frequency  branch  of  the  function  ft(u>,<p) 

(Eq.  4.52)  for  cp  »  60°.  For  a  given  u>,  the  slope 
of  the  curve  is  the  value  of  t/R  for  which  the 
phase  is  stationary  at  that  frequancy. 


fig.  5.2  -  Low  froquoncy  branch  of  tho  function  n(a>,(p) 

(lq.  4.52)  for  «p  »  60°.  For  a  g  ran  u>,  tho  slope 
of  tho  cur vo  la  tho  value  of  t/R  for  which  tho 
phoso  ia  stationary  at  that  froquoncy.  Roto  axis 
tonco  of  infloction  point,  at  which  tho  slopo 
has  its  minimum  volua. 
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approaches  8.  It  is  clear  that  the  function  has  an  inflection 

df} 

point  go  in  this  frequency  range.  The  derivative  decreases  from 

®  to  some  minimum  value  as  go  increases  from  8  cos  cp  to  go^  ,  and  then 

increases  again  to  ®  as  go  increases  from  go  to  8.  For  any  given  cp, 

cL  ^ 

the  minimum  value  of  represents  the  minimum  value  rf  t/R,  or  the 
maximum  value  of  R/t,  for  which  a  stationary  point  exists  (again 
only  for  the  exponential  with  the  minus  sign) .  For  values  of  R/t 
greater  than  the  maximum,  there  are  no  stationary  points  for  the 
given  cp,  whereas  for  values  smaller  than  the  maximum  there  are 
clearly  two  stationary  points.  At  the  maximum  value  of  R/t,  the 
two  stationary  points  coincide. 

This  situation  is  identical  to  that  which  arises  in  the  case  of 

incompressible  flow  in  a  density  stratified  liquid  considered  by 
(4) 

Mowbray  and  Rf.rity'  .  The  physical  consequences  are  the  same  as 
for  the  incompressible  case:  the  locus  in  physical  space  of  the 
double  points  of  stationary  phase  defines  a  front,  or  "caustic," 
representing  the  onset  of  tb.  disturbance.  Between  the  caustic  and 
the  spherical  front,  the  low  frequency  part  of  the  disturbance  is 
exponentially  small,  and  the  main  contribution  comes  from  the  high 
frequency  band. 

For  the  purposes  of  calculating  the  location  of  the  front,  it 
is  convenient  to  use  the  variables 


where  r  =  R  sin  cp  is  the  cylindrical  radius  and  z  =  R  cos  cp  is  the 
altitude.  In  terms  of  these  variables,  a  line  of  constant  phase 
5  ■  gu  -  —  n  can  be  written 

-flj  C*  .  .  fidi!  (5.4) 
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while  the  condition  for  stationary  phase  (5.2)  becomes 
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For  a  given  value  of  the  parameter  4,  which  corresponds  to  a 
particular  value  of  the  phase  at  a  given  time,  (5.4)  and  (5.5)  are 
parametric  equations  for  a  line  in  the  (£,T{) -plane  along  which  the 
particular  phase  is  stationary.  Each  point  on  such  a  line  corres¬ 
ponds  to  adifferent  value  of  ujf  the  point  of  stationary  phase.  A 
line  of  constant  and  stationary  phase  intersects  the  locus  of  double 
stationary  points  in  a  cusp.  The  condition  for  a  cusp  is 


du)  du) 


0  . 


(5.6) 


Applying  this  to  (5.5),  we  can  obtain  4  as  a  funct?  n  of  along  the 
front.  Equations  (5.4)  and  (5.5)  then  provide  a  one  to  one  mapping 
from  (|,u>)  to  (C,TD.  The  wave  front  obtained  in  this  manner  is  shown 
in  Fig.  5.3,  where  we  have  also  plotted  a  few  lines  of  constant  4 
for  the  purposes  of  illustration. 

The  co-ordinates  (C,TQ  of  any  point  on  the  front  are  the  hori¬ 
zontal  and  vertical  components,  respectively,  of  the  velocity  of  the 
front  at  that  point.  The  point  moving  along  the  ground  has  a  velocity 

“ft 

c/t  »  20  ■  ,9c  .  The  maximum  vertical  velocity  is  ,26c  ,  and  occurs 
at  an  angle  of  about  32°  with  the  horizontal.  (Every  point  on  the 
front  moves  radially  outward  from  the  origin.  ) 

As  just  indicated,  each  point  on  the  caustic  is  associated  with 
a  given  angle  cp  in  physical  space,  given  by  cp  =  tan  ^  T]/£.  Associated 
with  this  point  on  the  caustic  is  also  the  frequency  u>c  of  the  double 
point  of  stationary  phase.  This  is  the  first  frequency  (in  the  low 
frequency  band)  to  arrive  at  a  given  cp,  after  which  the  signal  splits 
into  two  frequency  components,  one  of  which  increases  with  time 
towards  0  and  the  other  of  which  decreases  towards  B  cos  cp.  In 
Fig.  5.4,  we  plot  in  polar  co-ordinates  the  frequency  as  a 

function  of  the  pole  angle  cp. 


30- 


r/t 


Fig.  5.3  -  Location  of  the  caustic  in  (r/t,  z/t)  space.  The 
dashed  lines  are  lines  of  stationary  phase  along 
vhich  the  phase  is  constant  at  any  instant. 
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rig.  5.4  -  Angular  frequency  at  tha  caustic,  in  unit* 
of  c*/h,  as  a  function  of  tha  pola  angla  <p. 
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In  light  of  the  above  discussion,  it  is  now  possible  to 
describe  briefly  the  frequency-time  history  of  the  signal  at  any 
location.  The  first  disturbance  to  arrive  at  any  location  (R,cp)  is 
the  acoustic  spherical  front,  containing  the  very  high  frequencies. 

As  time  progresses,  the  principal  frequency  in  this  part  of  the  dis¬ 
turbance  decreases  asymptotically  towards  Some  time  after 

the  passage  of  the  spherical  front,  the  caustic  arrives,  carrying 
the  frequency  associated  with  the  angle  cp  (  Fig.  5.4).  As  dis¬ 
cussed  above,  this  part  of  the  signal  then  splits  into  two  frequen¬ 
cies  ,  one  tending  toward  6  and  the  other  toward  0  cos  cp  as  time 
progresses.  Thus,  at  any  location,  the  signal  consists  asymptotically 
of  the  three  frequencies  8,  and  0  cos  cp.  This  information  is 
summarized  in  Fig.  5.5,  which  is  a  polar  graph  of  the  various 
groups.  Each  dashed  line  represents  the  location  in  (R/t,  0)  space 
of  a  group  of  a  given  frequency.  The  solid  line  is  the  caustic, 
which  appears  as  the  envelope  of  the  low  frequency  Dand.  Outside 
of  the  caustic,  only  the  high  frequency  band  is  present,  where# s 
inside  of  the  caustic  both  bands  clearly  contribute  to  the  signal. 


o 


Pig.  5. 


.£ 


5  -  Location  in  (R/t,cp)  space  of  groups  of  a  given 
frequency.  The  solid  l'ne  s  the  caustic,  which 
appears  as  the  envelope  of  the  low  freauencv  hand. 
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